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In order to shed light whether the 'even-odd conjecture' 
(even numbers of legs will superconduct accompanied by a 
spin gap while odd ones do not) for correlated electrons in lad- 
der systems, the pairing correlation is studied for the Hubbard 
model on a two- and three-leg ladders. We have employed 
both the weak-coupling renormalization group and the quan- 
tum Monte Carlo (QMC) method for strong interactions. For 
the two-leg Hubbard ladder, a systematic QMC (with a con- 
trolled level spacings) has detected an enhanced pairing corre- 
lation, which is consistent with the weak-coupling prediction. 
We also calculate the correlation functions in the three-leg 
Hubbard ladder and show that the weak-coupling study pre- 
dicts the dominant superconductivity, which refutes the naive 
even-odd conjecture. A crucial point is a spin gap for only 
some of the multiple spin modes is enough to make the ladder 
superconduct with a pairing symmetry (d-like here) compat- 
ible with the gapped mode. A QMC study for the three-leg 
ladder endorses the enhanced pairing correlation. 



I. INTRODUCTION 

Over the past several years, strongly correlated elec- 
tron systems with quasi-one-dimensional (ID) ladder 
structures have received much attention theoretically 
and experimentally. Experimental studies have received 
much impetus, since cuprate compounds-containing such 
structures have been fabricated recentlyE 

The_idea was inspired theoretically in 1986, when 
SchulzEl conjectured the following. If we consider a gas of 
repulsively interacting electrons on a ladder, the undoped 
system will be a Mott insulator, so that we may con- 
sider the system as an S — 1/2 antiferromagnetic (AF) 
Heisenberg magnet on a ladder. Then an AF ladder with 
A- legs should be similar to a AF £fj= N/2 single chain, 
which is exactly Haldane's systemJ3~l3 For the latter Hal- 
dane's conjecture predicts that the spin excitation should 
be gapless for a half-odd-integer spin (N: odd) or gap- 
ful for an integer spin (N: even). If the situation would 
be similar in ladders, a ladder having an even number of 
legs will have a spin gap, which should indicate that the 
ground state is a 'spin liquid' where the quantum fluctu- 
ation is so large that spins cannot order. On the other 
hand an odd number of legs will have gapless spin excita- 
tions, which should indicate that the ground state has an 
AF order. The presence of a spin gap in the former case 
may be a good news for superconductivity, since, there is 



a body of ideas dictating that a way to obtain supercon- 
ductivity is to carrier-dope a system that has a spin gap 
in the course of the study of high-Tc superconductivity. 
Such a scenario has been put forward by Rice et aln 

As far r, as. the spin gap is. .-concerned, both 
theoreticalurtLII and experimenta£3il3 studies on the 
undoped-ladder systems have indeed supported the con- 
jecture. 

The spin-gap-|Cpiijecture has recently been confirmed 
experimentallyllJlJ. Namely, a class of cuprates, 
Sr„_iCu„02n-i, has n-leg ladders on a Cu02 plane, 
and the two- leg ladders in SrCu203 exhibit a spin- liquid 
behavior characteristic of finite spin-correlation lengths, 
while the three-leg ladders in S^CuaOs have an AF be- 
havior. 

Rice et al. have further conjectured for doped systems 
that an even-numbered ladder should have a dominant 
interchain d-wave-like pairing correlation as expected 
from the persistent spin gap away from half-fillingJj The 
conjecture is partly based on an exact diagonalization 
study for finite systems for a two-leg t — J ladder by. 
Dagotto et aZ.^T-his-.was then followed by analyticaltj 
and numericalE30~Eil works on the doped t — J lad- 
der, which support the dominant pairing correlation in 
a certain region. In the phase diagram, the region for 
the dominant pairing correlation appears at lower val- 
ues of exchange coupling J than in the-, case of a sin- 
gle chain. Experimentally Uehara et alEU have recently 
observed superconductivity in a two-leg ladder material 
Sro.4Ca13.6Cu24O41.84 under high pressures. 

The Hubbard models on ladders are also of interest, 
since the Hubbard model may be regarded as an effec- 
tive model for cuprates (Fig.l). Since there is no exact 
solution for the Hubbard ladder, a most reliable anahd^. 
ical method at present is the weak-coupling theory£3'0 
which, in the continuum limit, linearizes the band struc- 
ture around the Fermi points to treat the interaction with 
a perturbative renormalization group. 




FIG. 1. Two- leg Hubbard ladder model; t(t±) and U are 
the intrachain(interchain) hopping and the on-site interac- 
tion, respectively. 



The weak-coupling theory has been applied to the two- 
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leg Hubbard ladder. BSEE At half-filling, the system 
reduces to-a spin-liquid insulator having both charge and 
spin gapsE3 with a finite SDW correlation length. Al- 
though one might expect that the Hubbard ladder would 
not exhibit a sizeable spin gap at half-filling unlike the 
t — J ladder, the spin gap for the Hubbard model es- 
timated with DMRG by Noack et alSMS is as large as 
Q.13t(~ 400K for t ~ 0.3eV for U = 8t with t± = t, 
which should correspond to the cuprates). The magni- 
tude of the spin gap is comparable with the spin gap (~ 
400K) experimentally estimated from the magnitude of 
susceptibility for SrCu203. 

When the carrier is doped, the weak-coupling theory 
supports the dominance of the pairing correlation whose 
symmetry is the same as that of the t — J ladder. The 
relevant scattering processes at the fixed point in the 
renormalization flow are the pair-tunneling process across 
the bonding and the anti-bonding bands (Fig. 2), and 
the backward-scattering process within each band. The 
importance of the pair-tunneling across the two bands 
(which exists in two or larger numbers of legs) for the 
dominance of pairing correlation in the two-leg Hubbard 
ladder is reminiscent of the Suhl-Kondo mechanism, that 
was proposed back in the 1950's for superconductivity, is, 
the transition metals with two (s- and d-like) bandsBJ'ta 
Muttalib and Emery have shown another example of 
the pair-tunneling mcchanism-for superconductivity with 
purely repulsive interactions!!!] More recently, the super- 
conductivity in t — t' — U model, which may be relevant 
for the chaipa that alternate with the ladder layers in 
the cupratea£j, has also been studied analytically^ and 
numericallyEa as a ID ladder-like system. 
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FIG. 2. Relevant pair-tunneling processes in two-leg Hub- 
bard ladder; Fig. (a) (Fig.(b)) is the forward (backward) type 
pair-tunneling process. 

The properties of the weak-coupling Hubbard ladder 
are similar to those of the t — J ladder for the regime 
where the pairing correlation is dominant: in addition 
to the existence of the spin gap, the duality relatioT&El, 
which suggests that the exponent for the pairing corre- 
lation (~ r 1 ^ 2Kp )) should be reciprocal to that for the 
Akp CDW correlation (~ r^ 2Kp )), holds in both .fchp-weak- 
coupling Hubbard ladder and the t— J ladder. EJ'ESEll Here 



K p is the critical exponent for the gapless charge mode, 
which tends to unity in the weak-coupling limit. The 
similarity may come from the form of the excitation gaps 
in the bosoniaation description. In both t — J and Hub- 
bard laddersEa, the only gapless mode is a charge mode 
with none of the spin modes being gapless ('OSO' phase 
in the language of the weak-coupling theorycll). 

However, there is a serious reservation for the weak- 
coupling theory. First, whether the model in the contin- 
uous limit is indeed equivalent to the lattice model is not 
obvious. More serious is the problem that the pertur- 
bational renormalization group is guaranteed to be valid 
only for an infinitcsimally small interaction strengths in 
principle. Specifically, when there is a gap in the exci- 
tation, the renormalization flows into a strong-coupling 
regime, so that the perturbation theory may well break 
down even for small interaction strengths. A way to check 
the reliability of the weak-coupling theory is to treat fi- 
nite systems with larger U with numerical calculations 
such as exact diagonalization, density-matrix renormal- 
ization gmw (DMRG), or quantum Monte Carlo (QMC) 
methodsE3'E3. Numerical calculations, on the other hand, 
have drawbacks due to finite-size effects. Thus the weak- 
coupling theory and the numerical methods should be 
considered as being complementary. 

Specifically, the dominance of the pairing correlation 
is indeed a subtle problem in numerical calculations. Ex- 
isting numerical resultsaSE^OO do appear to be con- 
troversial, where some of the results are inconsistent with 
the weak-coupling prediction as detailed in Sec. II. 

If we ignore these controversies, most of the exist- 
ing theories support the dominance of the pairing cor- 
relation in the doped two-leg ladders. Then, an even 
more important unresolved problem for superconductiv- 
ity in the doped ladders is the 'even-odd' conjecture. 
One can naively expect that the absence of spin gaps in 
odd-numbered legs will signify an absence of dominating 
pairing correlation, which will end up with an 'even-odd 
conjecture for superconductivity'. Indeed there has not 
been works looking into the pairing correlation functions 
for the three-leg ladder, which is the-simplest realization 
of odd- numbered legs. White et aZ.EJ have studied two 
holes doped in the t — J ladders at half-filling to find that 
two holes are bound in even-numbered (two or four) legs, 
while they are not in odd-numbered (three or five) legs, 
but (i)the existence of a binding energy for two carriers is 
not directly connected with the occurrence of supercon- 
ductivity, and (ii)the t-J model and the Hubbard model 
may exhibit different behaviors. 

Thus the second purpose of the present paper is to 
study the pairing correlation in Hubbard ladder models 
with three legs as compared with the two-leg case. 

The organization of the present paper is asJbllows. We 
first study the two-leg Hubbard ladder modelo having in- 
termediate interaction strengths with the QMC method 
in Sec. II. A new ingredient in this study is that we pay 
a special attention to the non-interacting (U = 0) single- 
particle energy levels which are discrete in finite systems. 
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We have found that if we make the levels on the two 
bands aligned (within an energy that is smaller than 
the spin gap) to mimic the thermodynamic limit, the 
weak-coupling result (an enhanced pairing correlation in 
the present case) is in fact reproduced for intermediate 
U. Further we check the effects of the inter- and intra- 
band Umklapp-scattering processes, which are expected 
to be present at-the special band fillings from the weak- 
coupling theoryEH. 

In Sec. Ill we study three(odd)-leg Hubbard ladder 
model (Fig. 3) with the weak-coupling thec*jfi2 using the 
enumeration of gapless modes by Arrigoni.E3 The system 
has one gapless spin mode and two gapful spin modes. 
Thus the gapful and the gapless spin modes coexist. The 
existence of the gapful modes is a result of the relevant 
interband pair-tunneling process across the top and the 
bottom bands (Fig. 4). We find that this gives rise to a 
dominant pairing correlation across the central and the 
edge chains reflecting the gapful spin modes, which coex- 
ists with the subdominant but power-law decaying SDW 
correlation reflecting the gapless spin mode. Thus the 
Suhl-Kondo-like mechanism for superconductivity exists 
not only in two-leg systems but also in a three(odd)-leg 
system as well, while, the SDW correlation also survives 
as expected. Schulzea also found similar results indepen- 
dently. 




FIG. 3. Three-leg Hubbard ladder model; t(t±) and U are 
the intrachain (interchain) hopping and the on-site interac- 
tion, respectively. 
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FIG. 4. Relevant pair-tunneling processes in the three-leg 
Hubbard ladder; Fig. (a) (Fig.(b)) is the forward (backward) 
type pair-tunneling process. 

We then study the three-leg Hubbard ladder model in 



Sec. IV with a QMC calculation^ as in the two-leg case. 
The technique to detect the enhanced pairing correlation 
in the two-leg case is also valid in the three-leg case. We 
found that the enhancement of the pairing correlation 
persists for the intermediate interaction strengths. We 
also study the effects of the Umklapp processes at special 
band fillings as in the two-leg case. 

The above numerical calculations for two or three legs 
are performed in a condition that the one-body energy 
levels are close to each other. We believe this condi- 
tion mimics the situation in the thermodynamic limit. 
In Sec.V, a circumstantial evidence for this is given by 
a QMC calculation of the pairing correlation in the ID 
attractive Hubbard model, which can be exactly solved 
with the Bethe Ansatzca, so that the exact asymptotic 
form of the correlation functions and the values of the 
spin gap are known at arbitrary values of parameters. 

II. QUANTUM MONTE CARLO STUDY OF THE 
PAIRING CORRELATION IN THE TWO-LEG 
HUBBARD LADDER 

In this section, we study the pairing correlation for the 
two-leg Hubbard ladder (Fig.l). The Hamiltonian of the 
two-leg Hubbard ladder is given in standard notations as 

H = -t^2(cl i(y c m+la + h.c.) 

aia 

-tx X^( c L'<r C2 ^°' + h -C.) + U ^n QiT n aii , (1) 

icr ai 

where a(= 1,2) specifies the chains and t(t±) is the 
intra(inter)-chain hopping. 

From an analytical point of view, if the system is free 
from Umklapp processes, the weak-coupling theory with 
the bosonization— fjombiuLcd with the renormalization- 
group techniquesc3'E3'L3E!l has indeed shown that the the 
two-leg Hubbard ladder has a spin gap and that the cor- 
relation function of the interchain d-wave-like pairing or- 
der parameter, d = {c\^C2ii — cuic 2 i^)/V^, decays as 
~ r~ 1 /( 2K >>} as a result of the relevant pair-tunneling pro- 
cess (Fig. 2), where K p is the critical exponent for the 
total-charge-density mode being only gapless and tends 
to unity in the weak-coupling limit. 

Since SDW and 2kp CDW correlations have to decay 
exponentially in the presence of a spin gap in a two- 
leg ladder, the only phase competing with pairing cor- 
relation will be 4&f CDW correlation, which should de- 
cay as r~ 2Kp . Hence the pairing correlation dominates 
over all the others if K p > 1/2. In numerical calcula- 
tions, however, the dominance of the pairing correlation 
in the Hubbard ladder appears to be a subtle problem. 
Namely, a DMRG study by Noack et al. for the doped 
Hubbard ladder with n — 0.875, U/t — 8, and t± = t 
shows no enhancement of the pairing correlation over the 
U = resaltcpS, while they do find an enhancement at 
t± = 1.5n2jO. Asai performed a quantum Monte Carlo 
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(QMC) calculation for-a 36 rungs ladder with n = 0.833, 
U/t = 2 and t± = 1.5tCj, in which no enhancement of the 
pairing correlation was found. On the other hand, Ya- 
maji et al. have found an enhancement for the values of 
the parameters when the lowest anti-bonding band lev- 
els for U = approach the highest occupied bonding 
band levels, although their results have not been-conclu- 
sive due to the small system sizes (< 6 rungs) o Thus, 
the existing analytical and numerical results appear to 
be controversial in the two-leg Hubbard ladder. 

Another point is that the above results are obtained 
away from special fillings where the Umklapp-scattermg 
processes are irrelevant. Recently, Balents and FisherEZl 
proposed a weak-coupling phase diagram (Fig. 5) which 
displays the numbers of the gapless spin and charge 
phases on the t± — n plane (where n is the band fill- 
ing) . The effects of the Umklapp processes are also dis- 
cussed there. At half-filling the interband Umklapp pro- 
cesses become relevant resulting in a spin-liquid insula- 
tor in which the pairing correlation decays exponentially. 
In addition, the intraband Umklapp process within the 
bonding band becomes relevant resulting in a gap in one 
charge mode in a certain parameter region where the 
bonding band is reduced to a half-filled band. This phase 
is called 'C1S2' phase because there are one gapless and 
two gapfull charge modes, while the phase at half-filling 
is called 'C0S0' phase because there is no gapless phase. 
We can expect that the pairing correlation decays ex- 
ponentially or is at least suppressed reflecting the exis- 
tence of the charge gap, although the direct calculation of 
the pairing correlation has not been done. Thus we also 
study the effects of the Umklapp processes in this section, 
keeping in mind the above weak-coupling results. 
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FIG. 5. Phase diagram in the weak-coupling limit (U — > 0) 
given in ref. [39]; the numbers of the gapless charge and spin 
modes (x and y, respectively) are denoted as CxSy and n is 
the band filling. In the dark region both of the two bands 
cross the Fermi level. 



the pairing in the Hubbard ladderEH with t± ~ t, es- 
pecially in order to clarify the origin of the discrepan- 
cies among the existing results. We conclude that the 
discreteness of energy levels in finite systems affects the 
pairing correlation enormously. 

The details of the QMC calculation are the follow- 
ing. We assumed the the periodic boundary conditions 
along the chain direction, cjv+i = c\ (where N labels 
the rungs) and took the non-interacting Fermi sea as the 
trial state The projection imaginary time r was taken to 
be ~ 60/t. We need such a large r to ensure the con- 
vergence of especially the long-range part of the pairing 
correlation. This sharply contrasts with the situation for 
single chains, where r ~ 20/t suffices for the same sam- 
ple length as considered here. The large value of r, along 
with a large on-site repulsion [/, makes the negative-sign 
problem serious, so that the calculation is feasible for 
U/t < 2. In the Trotter decomposition, the imaginary 
time increment [r/ (number of Trotter slices)] is taken 
to be < 0.1. We have concentrated on band fillings for 
which the closed-shell condition (no degeneracy in the 
non-interacting Fermi sea) is met. We set t = 1 in the 
remaining of this section. 

In the beginning we show in Fig. 6 the result for P(r) 
for t± = 0.98 and t± = 1.03 with U = 1 and the band 
filling n = 0.867 = 52 electrons/ (30 rungs x 2 legs). 
The U = result (dashed line) for these two values of t± 
are identical because the Fermi sea remains unchanged. 
However, if we turn on U , the 5% change in the tj_ = 
0.98 — > 1.03 is enough to cause a dramatic change in 
the pairing correlation: for t± = 0.98 the correlation 
has a large enhancement over the U = result at large 
distances, while the enhancement is not seen for t± = 
1.03. 
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In the remaining of this sectio n, ., w e perform an ex- 
tensive projector QMC calculationE3c3 to investigate the 
ground state correlation function P(r) = (Oj +r Oi) for 
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FIG. 6. Pairing correlation function, P(r), plotted against 
the real-space distance r in a 30 rungs Hubbard ladder having 
52 electrons for U = 1 with t± = 0.98 (□) and t± = 1.03 (o 
). The dashed line is the noninteracting result for the same 
system size, while the straight dotted line represents oc 1/r 2 . 
The solid line is a fit to the U = 1 result with t± = 0.98 (see 
text). 
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FIG. 7. Schematic image of the discrete energy levels of 
both bonding (0) and anti-bonding (n) bands for U — 0. 

In fact we have deliberately chosen these values to con- 
trol the alignment of the discrete energy levels at U — 0. 
Namely, when t± = 0.98, the single-electron energy lev- 
els of the bonding and anti-bonding bands for U = lie 
close to each other around the Fermi level with the level 
offset (Ae in Fig. 7) being as small as 0.004, while they 
are staggered for t± = 1.03 with the level offset of 0.1. 
On the other hand, the size .of. the spin gap is known to 
be around 0.05 for U — 8,oEj and is expected to be of 
the same order of magnitude or smaller for smaller values 
of U. The present result then suggests that if the level 
offset Ae is too large compared to the spin gap (which 
should be 0(0.01) for U ~ tp3i) the enhancement of the 
pairing correlation is smeared. By contrast, for a small 
enough Ae, by which an infinite system is mimicked, the 
enhancement is indeed detected in agreement with the 
weak-coupling theory, in which the spin gap is assumed 
to be infinitely large at the fixed point of the renormal- 
ization flow. Similar situation is also found in the case of 
a QMC-study of the pairing correlation in the t — t' — U 
model.E3 

Our-Jjesult is reminiscent of those obtained by Yamaji 
et al.,cB who found an enhancement of the pairing corre- 
lation in a restricted parameter regime where the lowest 
anti-bonding levels approach the highest occupied bond- 
ing levels. They conclude that the pairing correlation is 
dominant when the anti-bonding band 'slightly touches' 
the Fermi level. However, our result in Fig. 6 is obtained 
for the band filling for which no less than seven out of 
30 anti-bonding levels are occupied at U = 0. Hence the 
enhancement of the pairing correlation is seen to be not 
restricted to the situation where the anti-bonding band 
edge touches the Fermi level. 



However, one might consider that the enhancement of 
the pairing correlation in such a condition for the one- 
body energy levels should be rather due to finite size ef- 
fects, although we believe the condition is generally rel- 
evant for bulk systems. To clarify this point, we will 
further give a circumstantial evidence to justify that the 
condition for the one-body energy levels is relevant for 
bulk systems by a QMC study for ID attractive Hub- 
bard model in Sec.V. 

Now, let us more closely look into the form of P(r) 
for t± = 0.98. It is difficult to determine the exponent 
from results for finite systems, but here we attempt to fit 
the data by assuming a trial function expected from the 
weak-coupling theory. Namely, we have fitted the data 
with the form, 

p W = iEK 1/2 

d=± 

+ [(2 - c) - cos(2fc° r d ) - cos(2^r d )]r- 2 } (2) 

with the least-square fit (by taking logarithm of the 
data). Because of the periodic boundary condition, we 
have to consider contributions from both ways around, 
so there are two distances between the 0-th and the r-th 
rung, i.e, r + — r and r_ = N — r. The periods of the 
cosine terms are assumed to be the non-interacting Fermi 
wave numbers of the bonding and the anti-bonding bands 
in analogy with the single-chain case. 

The overall decay should be 1/r 2 as in the pure ID case 
in the weak-coupling limit. We have assumed the form 
c/r 1 / 2 as the dominant part of the correlation at large 
distances because this is what is expected in the weak- 
coupling theory. Here c is the only fitting parameter in 
the above trial function. A finite U ~ 1 may give some 
correction, but the result (solid line in Fig. 6) fits to the 
numerical result surprisingly accurately with a best-fit 
c = 0.10. If we least-square fit the exponent itself as 
l/r Q , we have a < 0.7 with a similar accuracy. Thus 
a finite U may change a, but a > 1 may be excluded. 
To fit the short-range part of the data, a non-oscillating 
(2 — c)/r 2 term is required, which is not present in the 
weak-coupling theory. We believe that this is because the 
weak-coupling theory only concerns with the asymptotic 
form of the correlation functions. 

In Fig. 7, we show a result for a larger system size (42 
rungs) for a slightly different electron density, n = 0.905 
with 76 electrons and t± = 0.99. We have again an ex- 
cellent fit with c = 0.07 this time. 
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FIG. 8. Plot similar to that in Fig. 6 for a 42 rungs system 
having 76 electrons with t± = 0.99. 



In Fig. 8, we display the result for a larger U = 2. 
We again have a long-ranged P(r) at large distances, 
although P{r) is slightly reduced from the result for 
[7=1. This is consistent with the weak-coupling the- 
ory, in which K p is a decreasing function of U so that 
after the spin gap opens for U > 0, the pairing correla- 
tion decays faster for larger values of U. 




FIG. 9. Plot similar to that in Fig.6 for U = 2. 



Now we explore the effects of the Umklapp processes. 
For that purpose we concentrate on the filling depen- 
dence for a fixed interaction U = 2. We have tuned the 
value of t± to ensure that the level offset (Ae) at the 
Fermi level is as small as 0(0.01) for U = 0. In this way, 
we can single out the effects of the Umklapp processes 
from those due to large values of Ae. If we first look at 
the half-filling (Fig. 10), the decaying form is essentially 



similar to the U = result. At half-filling, the inter- 
band Umklapp processes emerge and, according to the 
weak-coupling theory, open a charge gap, which results 
in an exponential decay of the pairing correlation. (We 
should note that there are two kinds of charge gaps. The 
one, which is produced by the pair-tunneling processes, 
causes the long range order of the Josephson phase result- 
ing in the enhancement of the pairing correlation, while 
the other, which is produced by the Umklapp processes, 
causes the long range order of the phase of the CDW 
resulting in the suppression of the pairing correlation) 




FIG. 10. Pairing correlation P(r) (□) against r for a 30 
rungs system for U = 2 with t± = 0.99 and 60 electrons 
(half-filled). The dashed line represents the non-interacting 
result. 
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FIG. 11. Pairing correlation P(r) (□) against r for a 30 
rungs system for U = 2 with t± = 1.01 and 40 electrons 
(half- filled bonding band). The dashed line represents the 
non-interacting result. 
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It is difficult to tell from our data whether P(r) decays 
exponentially. This is probably due to the smallness of 
the-jdiarge gap. In fact, the DMRG study by Noack et 
aZ.EjH have detected an exponential decay for larger val- 
ues of U, for which a larger charge gap is expected. 

When n is decreased down to 0.667 (Fig. 11), we again 
observe an absence of enhancement in P(f\- This is again 
consistent with the weak-coupling theorycJ: for this band 
filling, the number of electrons in the bonding band co- 
incides with N{= 30) at U = 0, i.e., the bonding band 
is half-filled. This will then give rise to intraband Umk- 
lapp processes within the bonding band resulting in the 
'C1S2' phase as discussed in Section 2.1. The spin gap 
is destroyed and the singlet-pair of electrons or holes are 
prevented from forming, so that the pairing correlation 
will no longer decay slowly there. Noack et aVEi have 
suggested that the suppression of the spin gap and t 
pairing correlation function around t± = OAt in ref. 
may be due to the intraband Umklapp process. 

In this section, we have detected the enhancement of 
the pairing correlation which is consistent with the weak- 
coupling theory. 

We have also seen that there are three possible causes 
that reduce the pairing correlation function in the Hub- 
bard ladder: 

(i) the discreteness of the energy levels, 

(ii) reduction of K p for large values of U/t, and 

(iii) effect of intra- and interband Umklapp processes 
around specific band fillings. 

The discreteness of the energy levels is a finite-size ef- 
fect, while the others are present in infinite systems as 
well. We can make a possible interpretation for the exist- 
ing results in terms of these effects. For 60 electrons on 
36 rungs with t± — 1.5t in refo, for instance, the non- 
interacting energy levels have a significant offset ~ 0.15i 
between bonding and anti-bonding levels at the Fermi 
level, which may be the reason why the pairing correla- 
tion is upt enhanced for U/t — 2. For a large U/t(= 8) in 
ref£3o€3, (ii) and/or (iii) in the above may possibly be 
important in making the pairing correlation for t± = t 
not enhanced. The effect (iii) should be more serious for 
t± = t than for t± = l.5t because the bonding band is 
closer to the half-filling in the former. On the other hand, 
the discreteness of the energy levels might exert some ef- 
fects as well, since the non-interacting energy levels for a 
32-rung ladder with 56 electrons (n = 0.875) in an open 
boundary condition have an offset of 0.15t at the Fermi 
level for t± = t while the offset is 0.03i for t± = 1.5t. 

Let us comment on a possible relevance of the present 
result to the superconductivity reported recently for a 
cuprate laddered, especially for the pressure dependence. 
The material is Sro.4Ca13.6Cu24O41.s4, which contains 
layers consisting of two-leg ladders and those consisting 
of ID chains. Superconductivity is not observed in the 
ambient pressure, while it appears with Tq ~ WK under 
the pressure of 3 GPa or 4.5 GPa, and finally disappears 
at a higher pressure of 6 GPa. This material is doped 
with holes with the total doping level of 6 = 0.25, where 



5 is defined as the deviation of the density of electrons 
from the half-filling. It has been proposed that at ambi- 
ent pressure the holes are mostly in the chains, while high 
pressures cause the carrier to transfer into the laddersc3. 
If this is the case, and if most of the holes are trans- 
ferred to the ladders at 6 GPa, the experimental result 
is consistent with the present picture, since there is no 
enhancement of the pairing correlation for 5 = and 

6 ~ 0.3 due to the Umklapp processes as we have seen. 
Evidently, further investigation especially in the large-?/ 
regime is needed to justify this speculation. 



III. WEAK-COUPLING STUDY OF 
CORRELATION FUNCTIONS IN THE 
THREE-LEG HUBBARD LADDER 

In this section, we study correlation functions for the 
three(odd)-leg Hubbard ladder by the weak-coupling the- 
ory. As discussed in Sec. I, an increasing fascination in 
ladder systems has been caused by an 'even-odd' conjec- 
ture for the existence_of spin gap by Schulztj and indepen- 
dently by Rice et alu at half-filling. When the system is 
doped with carriers, it is naively supposed that an even- 
numbered ladder should exhibit superconductivity with 
the interchain singlet pairing as expected from the per- 
sistent spin gap, while an odd-numbered ladder should 
have the usual 2k p SDW reflecting the gapless spin ex- 
citations. 

Theoretically, however, whether the 'even-odd' conjec- 
ture for superconductivity continues to be valid for triple 
chains remains an open question. There had been no re- 
sults for the pairing correlation function in the three-leg 
t - J or Hubbard ladder (Fig.3). 

On the other hand, Arrigoni has looked into a three- 
leg with weak Hubbard-type interactions by the usual 
perturbational renormalization-group technique, which is 
quite similar to that developed by Balents and Fisher for 
the two-leg case,EZl to conclude that gapless and gapful 
spin excitations coexist there.E-3 

Namely, he has actually enumerated the numbers of 
gapless charge and spin modes on the phase diagram 
spanned by the doping level and the interchain hopping. 
He found that, at half- filling, one gapless spin mode ex- 
ists for the interchain hopping comparable with the in- 
trachain hopping, in agreement with some experimental 
results and theoretical expectations (Fig. 12). Away from 
the half-filling, on the other hand, one gapless spin mode 
is found to remain at the fixed point in the region where 
the fermi level intersects all the three bands in the non- 
interacting case. From this, Arrigoni argues that the 2k p 
SDW correlation should decay as a power law as expected 
from experiments. 
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FIG. 12. Phase diagram in the weak-coupling limit 
(U — > 0) given in ref. [53]; the numbers of the gapless charge 
and spin modes (x and y, respectively) are denoted as CxSy 
and p = 2 — 7i (n is the band filling). In the dark region all 
the three bands cross the Fermi level. 



On the other hand, his result also indicates that two 
gapful spin modes exist in addition. While a spin gap 
certainly favors a singlet superconducting (SS) correla- 
tion when there is only one spin mode, we are in fact 
faced here with an intriguing problem of what happens 
when gapless and gapful spin modes coexist, since it may 
well be possible that the presence of gap(s) in some out 
of multiple spin modes may be sufficient for the domi- 
nance of a pairing correlation. Furthermore, as discussed 
below, the gaps of two spin modes emerge as an effect of 
the pair-tunneling process across the top and the bottom 
bands (Fig. 4). This is remiruscemt-pf the two-leg case and 
of the Suhl-Kondo mechanism. These have motivated 
us, in this section, to actually look at .the correlation func- 
tions using the bosonization methods at the fixed point 
away from half-filling. Although in the three-leg case, 
we can consider the two boundary conditions across the 
legs, i.e., open boundary condition (OBC) and periodic 
boundary condition (PBC), here we concentrate on the 
open boundary condition (OBC) across the chains, where 
the central chain is inequivalent to the two edge chains. 
The reason is that we would like to (i) study the realis- 
tic boundary condition corresponds to cuprates, and (ii) 
to avoid the frustration introduced in the periodic three- 
legs. 

We find that the interchain SS pairing across the cen- 
tral and edge chains is the dominant correlation, which 
is indeed realized due to the presence of the two gap- 
ful spin modes. On the other hand, the SDW corre- 
lation, which has a slowly-decaying power law for the 
intra-edge chain reflecting the .gapless spin mode.-coex- 
ists but is only subdominant.a Recently SchulzCJ has 
independently shown similar results for a subdominant 
2Icf SDW and the interchain pairing correlations which 
are given in this section. 



A. Model and the Calculation 

The three-leg Hubbard model with OBC is defined by 
the following Hamiltonian, 

H = -tY^( c l ia c m+i° + h - c -) 

~t± 5^(4ia c /K" + c\ ia C lw + h.C.) 

+ u ^2 n ^l n n*l' ( 3 ) 

where t(t±) is the intra-(inter-)chain hopping, i labels the 
rung while fi = a, /?, 7 labels the leg (with (3 being the 
central one). In the momentum space we have 



H = ^2 f-2*cos(fc) - V2i_L j a\ ka a lkl7 

ka 

-2t 2J cos{k)a\ kl7 a 2 ka 

ka 

+ (-2icos(fc) + V2tj_) al k(j a 3k o 

ka 

-\-U (terms of the form a^a^aa). 



(4) 



Here ajka annihilates an electron with lattice momentum 
k in the j-th band (j — 1,2,3), where a^ka is related 
to c^ka (the Fourier transform of c^a) through a linear 
transformation, 







Cf3ka 


"I 












1 a 2 ka 


(5) 


) V a 3ka J 





Hereafter we linearize the band structure around the 
fermi points as usual and neglect the difference in the 
fermi velocities of three bands, as is done for calculating 
the correlation functions di re c tly in the weak-coupling 
theory for the two-leg case,c§tlil which will be accept- 
able for the weak interchain hopping. These approxi- 
mations enable us to calculate the correlation functions. 
The difference in fermi velocities of three bands will not 
be important qualitatively as long as we consider the 
case where three bands cross the fermi energy, for which 
Arrigoni's result falls on the same strong coupling fixed 
point on the plane of interchain hopping and filling. In 
the following, we focus on the case in which all of three 
bands are away from half-filling. 

The part of the Hamiltonian, H^, that can be diagonal- 
ized in the bosonization only includes forward-scattering 
processes in the band picture, and has the form 

-^d = H S p[ n -)- H charge , 

ff S pi„ = J2^f dx[j^(d x <t> i+ ) 2 + K m (d x &.) 2 ], (6) 
u r \ 
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Here </>j+ is the spin phase held of the i-th band, Xi+ 1S 
the diagonal charge phase held, while (xi- ) is the held 
dual to 0i+(xi+)i Ktri{Kpi) the correlation exponent for 
the 4>{xi) phase with v„i{v P i) being their velocities. For 
the Hubbard-type interaction, we have v a i — vf, K„i=l 
for all i's, while v p \ = vf, v p i — «f\A — 4g 2 , v p3 — 
v F yi- g 2 /^ K pl = 1, K p2 = - 2g)/(l + 2g), 

K p3 = v /(l- 5 /2)/(l + ff /2), where g = U/2irv F is the 
dimensionless coupling constant. The derivation of the 
above equation is given in the Appendix. 

The diagonalized charge held Xi± 1S linearly related to 
the initial charge held 8i± of the i-th band as 







(a) 


■( 














(7) 







where both 9i± and <f>± are related to the held operator 
for electrons i/ , i+(-)o'j which annihilates an electron on 
the right- (left-) going branch in band i as 

4>i+(-)<?( x ) = e^p{±ik iF x 

± l -[6 l+ {x) ± 6i_(x) + a(4> i+ (x) ± &_(a:))]}. 

Here the r\i ra 's are Haldane's U operators^ which ensure 
the anti-commutation relations between electron opera- 
tors through the relation, {r)i ra -, rii'r'a'}+ = ^n'&rr'baa 1 •, 

Vira = Vira- 

There are still many scattering processes correspond- 
ing to both the backward and the pair-tunneling scat- 
tering processes, which cannot be treated exactly. Ar- 
rigoni examined the effect of such scattering processes 
by the perturbational renormalization-group technique. 
He found that the backward-scattering interaction within 
the hrst or the third band turn from positive to negative 
as the renormalization is performed and that the pair- 
tunneling processes across the hrst and third bands also 
become relevant. As far as the relevant scattering pro- 
cesses are concerned, the hrst (third) band plays the role 
of the bonding (anti-bonding) band in the two-leg case. 
At the hxed point the Hamiltonian density, H* , then 
takes the form, in term of the phase variables, 



H* 



- |gl C os(20 1+ (x)) - Mgco S (20 3+ (z)) 
2g/t(1 ' 3) cos(V2 Xl _(z))sin0 1+ (z)sin0 3+ (x) , (9) 



+ 



7T 2 A 2 



where both <7f,(l) and c/f,(3) are negative large quantities, 
and gf t (l,3) is a positive large quantity. 

This indicates that the phase fields 4>\+, <f>3+, and Xi- 
are long-range ordered and hxed at tt/2, tt/2, and tt/V2, 
respectively, which in turn implies that the correlation 
functions that contain </>3_, and Xi+ helds decay 

exponentially. The renormalization procedure will affect 
the velocities and the critical exponents for the gapless 



helds, X2±: X3±? and 4>2±, so that we should end up with 
renormalized w*'s and if*'s. 

In principle, the numerical values of renormalized u*'s 
and K*'s for hnite g may be obtained from the renor- 
malization equations as has .been attempted for a double 
chain by Balents and FisherEH, although it would be dif- 
ficult in practice. However, at least in the weak-coupling 
limit, g — > 0, to which our treatment is meant to fall 
upon, we shall certainly have v* ~ vf and K* ~ 1 for 
gapless modes even after the renormalization procedure. 



B. Results for the Correlation Functions 

Now we are in position to calculate the correlation 
functions, since the gapless helds have already been di- 
agonalized, while the remaining gapful helds have the 
respective expectation values. The details of the calcu- 
lation of the correlation functions are given in the Ap- 
pendix. The two-particle correlation functions which in- 
clude the following two particle operators in the band 
description are shown to have a power-law decay: 

(1) operators constructed from two operators involving 
only the second band (since the charge and the spin 
!p)hases are both gapless, electrons in this band should 
have the usual TL- liquid behavior), 

(2) order parameters of singlet superconductivity within 
the hrst or third band, Vi+tU)^i-i(T)' ^3+TU)^3-|(T)- 
As a result, the order parameters that possess power-law 
decays should be the following, 

(A) The correlations within each of the two edge (a and 
7) chains or across the two edge chains: 

(a) 2k F CDW, 

Ointra2fc F CDW = V'l( 7 ) + | ^a( 7 )-T ! 
t 

a(7)4 

(b) 2k F SDW, 



OintcrCDW — ^l(-,) + ^l(a)-^ , 



intraSDW 



^a( 7 )+T^«(T)-l5 



O 

OintcrSDW = ^i( 7 ) + |'07(a)-J.) 

(c) singlet pairing (SS), 

OintraSS = ^ a^ + ^a^)- 
OintorSS = '/ , a(7) + T'^7(a)-l' 

(d) triplet pairing (TS), 

OmtraTS = 1pa(^) + ]i>a(^)-] ') 
OintcrTS = ^a(^) + ] ^^(a)-] > 

(B) The 4fci? CDW which is written with four electron 
operators, 

O^cdw = ipl +1 ipl +L tpu~^tp^-i (v = a, (3, 7), 

(C) The singlet pairing across the central chain (/3) and 
edge chains (Fig. 13), 

OCESS = E CT cr (V'a+<T + Ipj+a^p-.-a- 







a 
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FIG. 13. Schematic picture for the interchain (CESS) pair- 
ing in the doped three-leg Hubbard ladder. 

In the band picture we can rewright Ocess as compris- 
ing 



OCESS ~ <7(i)l+<ri>l 



(10) 



We cannot easily name the symmetry of the pairing, 
although we naively might call this pairing d-wave-like 
in a similar sense as in the two-leg case, in which a 
pair is called d-wave when the pairing, in addition to 
being off-site, consists of a bonding b and .aa d . an anti- 
bonding band pairs with opposite signs .EZH23'E3'E2J , E2i Thus 
the edge-chain SDW correlation has a power-law decay, 
while the SDW correlation within the central chain de- 
cays exponentially since it consists of the terms contain- 
ing 4>x- and/or </>3_ phases. Although we consider the 
case away from half-filling, the SDW correlation should 
obviously be more enhanced at half-filling. The experi- 
ments at half-filling do not contradict the present results, 
since the experiments should detect the total SDW corre- 
lation of all the chains and the SDW correlation is more 
enhanced at half-filling. However the present theory cor- 
responds only to the infinitcsimally small interaction in 
principle, although the actual cuprates have a strong in- 
teractions between electrons. 

Intra- or inter-edge correlation functions have to in- 
volve forms bilinear in a-ika in eq.(|^). They are described 
in terms of the charge field 82 for the second band, which 
does not contain xi, a phase-fixed field (see eq.(0)). Thus 
the edge-channel correlations are completely determined 
by the character of the second band (the Luttinger-liquid 
band), while the other phase fields, being gapful, are irrel- 
evant. The final result for the edge-channel correlations 
at large distances, up to 2kp oscillations, is as follows re- 
gardless of whether the correlation is intra- or inter-edge: 

(O 2feF CDwWOL F cDw(0)> ~ x~^ K '^U)-\ 
(Osdw(x)oLw(O)) ~ x-^ +2K U)'\ 

(OssWOL(0))~^ §( % +7 ^ ) " 1 , (11) 

(where we have put i-Cl = 1 for the present spin- 
independent interaction^) In addition, the 41cf CDW 
correlation decays as 



By contrast, if we look at the pairing Ocess {%) across 
the central chain and the edge chains, this pairing, which 
circumvents the on-site repulsion and is linked by the res- 
onating valence bonding across the neighboring chains, is 
expected to be stronger than other correlations as in the 
two-leg case. The correlation function for Ocess ( x ) is 
indeed calculated to be 



(Ocess (z)O cess (0)> 



3 I k' + : 
X P 2 



-) 



(13) 



From the calculations given in the Appendix, we can see 
that the interchain pairing exploits the charge gap and 
the spin gaps to reduce the exponent of the correlation 
function, in contrast to the intra-leg pairing. In addition 
to that, we also find that the roles of the first (third) 
band corresponds to those of the bonding (anti-bonding) 
band in the two-leg CcLSG, clS far as the dominant pair- 
ing correlation is concerned. If we consider the weak- 
interaction limit (U — > +0) as in the two-leg case, all 
the K*'s will tend to unity, where the CESS correlation 
decays as a; -1 / 2 while those of other correlations decays 
as x~ 2 at long distances. Thus, at least in this limit, the 
CESS correlation dominates over the others. The duality 
(which dictates that the pairing and density-wave expo- 
nents are reciprocal of each otheiElJ) is similar to that in 
the two-chain case, in which the interchain SS decays as 
x ' while that of the Akp CDW decays as x 

In this section, we have studied correlation functions 
using the bosonization method at the renormalization- 
group fixed point, which was obtained by Arrigoni, away 
from half-filling in the region where the fermi level in- 
tersects all the three bands in the non-interacting case. 
We found that the interchain singlet pairing across the 
central chain and either of the edge chains is the domi- 
nant correlation contrary to the naive 'even-odd' conjec- 
ture for the superconductivity in ladder systems, while 
the SDW correlations in two edge chains coexist but are 
subdominant. The power law decay of the SDW corre- 
lation does not contradict with the even-odd conjecture 
at the half-filling, where the Umklapp scattering play a 
important role resulting in an enhancement of the SDW 
correlation. 

The renormalization study is valid only for in- 
finitesimally small interaction strengths and sufficiently 
small interchain hoppings in principle, while the actual 
cuprates have strong interactions between electrons, so 
that the relevance of the present results to the real ma- 
terials is uncertain. However the present study suggests 
an important theoretical message that the dominance of 
superconductivity only requires the existence of gap(s) 
in some spin modes when there are multiple modes in 
multi-leg ladder systems no matter whether the number 
of legs is odd or even. 



(O 4fcFC Dw(x)ol fcFCDW (0)) ~ x -i(^; 2+ K; 3 )_ (12) 
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IV. QMC STUDY OF THE PAIRING 
CORRELATION IN THE THREE-LEG 
HUBBARD LADDER 

In the previous section, we have discussed the correla- 
tion functions in the three-leg Hubbard ladder within the 
weak-coupling theory. A key point in the previous sec- 
tion is that the gapless and the gapful spin excitations 
coexist in a three-leg ladder and the modes give rise to a 
peculiar situation where a specific singlet pairing across 
the central and edge chains (CESS pairing), that may 
be roughly a d-wave pairing, is dominant, while the 2kp 
SDW on the edge chains simultaneously shows a subdom- 
inant but still long-tailed (power-law) decay associated 
with the gapless spin mode. This result is stimulating 
since it serves as a counter-example of a naive 'even-odd' 
conjecture. 

However, there is a serious question about these weak- 
coupling results as discussed in the two-leg case in Sec. II. 
First, only for infinitesimally small interactions and suf- 
ficiently small hoppings are the results in Sec. Ill guaran- 
teed to be valid in principle. Furthermore, when there is 
a gap in the excitation, the renormalization flows into a 
strong-coupling regime, so that the weak-coupling theory 
might break down even for small U . Hence it is impera- 
tive to study the problem from an independent numeri- 
cal method for an intermediate strength of the Hubbard 
U ~ t and an interchain hopping t± ~ t. Although such 
a comparison of the numerical result for U *~ t with the 
weak-coupling theory has been done for the two-leg sys- 
tem in Sec. II, this does not necessarily shed light on the 
situation in the three-leg case, where gapless and gapful 
modes coexist. In this section, it is shown that the QMC 
result for the three-leg Hubbard ladder indeed turns out 
to exhibit an enhancement of the pairing coo-elation even 
for finite coupling constants, U/t = 1 ~ 2.c3 

In addition, we also study the effects of various Umk- 
lapp processes at special fillings as in the|-t,wo-leg case 
keeping in mind the above Arrigoni's worked which also 
studied the effects of some Umklapp processes with the 
weak-coupling theory. 

Throughout this section, we concentrate on the case in 
which all three bands cross the Fermi surface to explore 
the properties of a three-band system. .— . .— . 

The projector Monte Carlo method is employec£j'E£l 
to investigate the ground-state pairing correlation func- 
tion P(r) = {0]O j+r ), where O t = C ESs{i) = (c aia + 
c-yio-)cpi-cr— (e m -_ CT + e 7 i_ a -)c / a, a -. We assume the periodic 
boundary condition along the chain direction, cjv+i = ci, 
where N is the number of rungs. We first consider the 
case where the intra- and the inter-band Umklapp pro- 
cesses are irrelevant because that is where the above- 
mentioned weak-coupling theory is valid. The details of 
the QMC calculation are similar to those for our QMC 
study for the two-leg case. Specifically, the negative-sign 
problem makes the QMC calculation feasible for U <2t 
as in the two-leg case. We set t = 1 hereafter. 



In the two- leg case with a finite U, we have found an 
interesting property for finite systems: the pairing corre- 
lation is enhanced in agreement with the weak-coupling 
theory only when the single-electron energy levels of the 
bonding and the anti-bonding bands lie close to each 
other around the Fermi level (which is certainly the case 
with an infinite system). When the levels are misaligned 
(for which a 5% change in t± is enough), the enhance- 
ment of the pairing correlation dramatically vanishes. In 
the weak-coupling theory, the ratio of the spin gap to 
the level offset is assumed to be infinitely large at the 
fixed point of the renormalization flow, so that the spin 
gap should naturally be detectable in finite systems only 
when the level offset is smaller than the gap. 

We have found that this applies to the three-leg ladder 
as well, i.e., the pairing correlation is enhanced when 
the single-electron levels of the first and third bands lie 
close to each other. Hence we concentrate on such cases 
hereafter. 

In the beginning we show in Fig. 14 the result for P(r) 
for t± = 0.92 with U = 1 with the band filling n = 
0.843 = 86 electrons/(34 rungsx 3 sites). For this choice 
of t± the levels in the first and the third bands lie close 
to each other around the Fermi level within 0.01. We 
can see that a large enhancement over the U — result 
at large distances indeed exists. This is the key result of 
this section. 




FIG. 14. QMC result for the pairing correlation func- 
tion, P(r)(D), plotted against the real space distance r in 
a three-leg Hubbard ladder with 34 rungs having 86 electrons 
for U — 1 with t± = 0.92. The dashed line is the noninteract- 
ing result for the same system size, while the straight dashed 
line represents ~ r~ 2 . The solid line is a fit to a trial function 
(see text). 

Although it is difficult to determine the decay exponent 
of P{r), we can fit the data by supposing a trial function 
as expected from the weak-coupling theory as we did in 
the two-leg case, 
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d=± 
+[(2-c)- 



1/2 



cos(2fcFir d ) - cos(2k F3 r d )}r d 2 } (14) 



Here £;_Fi(fcF3) is the non-interacting Fermi wave num- 
ber of the first (third) band, while a constant c, which 
should vanish for U = 0, is here least-square fit (by tak- 
ing logarithm of the data) as c = 0.05. As in the two-leg 
case, since we assume the periodic boundary condition, 
we have to consider contributions from both ways around, 
so there are two distances between the 0-th and the r-th 
rung, i.e., r + = r and r_ = N — r. The overall decay 
should be 1/r 2 as in the single-chain case, while the term 
c/r 1 / 2 , £ ne dominant correlation at lar ge ,, d istances, is 
borrowed from the weak-coupling result .c§c3 The QMC 
result for a finite U = 1 fits to the trial form (solid line in 
Fig. 14) surprisingly accurately. A finite U may give some 
corrections to these functional forms, but even when we 
best- fit the exponent itself as c/r a in place of c/r 1 / 2 , we 
obtain a < 0.7 with a similar accuracy. 

In Fig. 15, we show the result for a larger interaction 
U = 2. The result again shows an enhanced pairing 
correlation at large distances. However, the enhancement 
is slightly reduced than that in the U — 1 case. This is 
consistent with the weak-coupling theory, in which K* 's 
should decrease with U. 




FIG. 15. Plot similar to that in Fig. 14 for U = 2. 



10" 



10" 



10" 



fl\M 



band 1 2 3 



t| ,o.oi v V \/ V 



10 u 



10' 



FIG. 16. Similar plot as in Fig. 14 for a 38 rungs system 
having 82 electrons for U = 1 with t±_ — 0.685. The inset 
schematically depicts the positions of energy levels for the 
noninteracting case. 



Furthermore, we study if the presence of the second 
band around Ep can be detrimental to superconductiv- 
ity. In Fig. 16, we make the single-electron energy levels 
of all the three bands lie close to each other around the 
Fermi level. This is accomplished here for t± = 0.685 
and the band filling n = 0.719 = 82 electrons/(38 rungs 
x 3 sites). The highest occupied level of the second band 
then lies between that of the first band and the lowest 
unoccupied level of the third band (lying above the high- 
est occupied level of the first band by as small as 0.01, 
inset of Fig. 16). 

The result in Fig. 16 for U — 1 shows that the pairing 
correlation is enhanced as well. Thus we may consider 
that the second band does not hinder the enhancement 
of the pairing correlation in other bands. This is also 
consistent with the weak-coupling theory, in which all of 
the scattering processes connected with the second band 
are irrelevant. The fit of the correlation function to the 
trial one is again excellent with c = 0.03. 

In the remaining of this section, we discuss the effects 
of the Umklapp processes at special fillings to clarify the 
doping dependence. In the three-leg Hubbard model, the 
Umklapp processes can play an important role at specific 
band fillings. 

ArrigonC-l also studied the effect of Umklapp processes 
within the weak-coupling theory, although he did not cal- 
culate the correlation functions directly (see Fig. 12). He 
studied two cases that have the relevant Umklapp pro- 
cesses. 

(i) the half-filled case. 

(ii) the case when the bottom band is half-filled, in 
which the intraband Umklapp process may become rele- 
vant within the first band. 

In both cases, the Umklapp processes become relevant. 
In the former case, the system has a gapless spin exci- 
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tation suggesting a power-law decaying AF correlation, 
as discussed by Arrigoni, for t± <~ t, a region of inter- 
est. The phase is called a 'C1S1' regime, since there is 
one gapless mode in charge or spin mode. Although the 
existence of a gapless charge mode suggests that the sys- 
tem is not an insulator, a full charge gap is expected to 
appear for sufficiently large U. Then, the pairing corre- 
lation will be suppressed, although the direct calculation 
has not been done. (As discussed in Sec. II, there are two 
kinds of charge gaps. The one, which is produced by 
the pair-tunneling processes, favors the pairing correla- 
tion as in the previous section, while the other, which is 
produced by the Umklapp processes, suppress the pairing 
correlation as in the present section.) 

In the latter case, the Umklapp process also becomes 
relevant resulting in a 'C2S3' phase with the two gap- 
less charge modes and three gapless spin modes and the 
singlet-pair of holes or electrons is prevented from bind- 
ing. Thus the pairing correlation will be not enhanced 
because of the absence of gapful spin mode(s). This sit- 
uation is reminiscent of the two-leg case, in which the 
spin gap is destroyed when the bonding band is half- 
filled. The Umklapp process within the first band in the 
present case is identified with that in the bonding band in 
the two-leg case, since the first (third) band corresponds 
to the bonding (anti-bonding) band as far as the pairing 
correlation is concerned as discussed in Sec. III. 

Given this situation, our motivation here is to look 
at the pairing correlation and, in addition, to explore 
in which case the Umklapp scattering does or does not 
affect the pairing correlation. We study the above two 
cases and also study the case in which the second band 
is half-filled. Namely, we wish to see whether both the 
first and the third bands, which involve the pairing or- 
der parameter, are affected by an indirect effect of the 
Umklapp process in the second band. Such a situation 
emerges in ladder systems with three or larger number of 
legs. 

We have tuned the value of t± to ensure that the level 
offset (Ae) between the first and the third bands at the 
Fermi level is as small as 0(0.01) for U = to single 
out the effect of Umklapp processes from those due to 
large values of Ae. Ideally, we should make the single 
electron energy levels of all the three bands lie close to 
each other around the Fermi level but that is impossible 
within the tractable system sizes. However, when the 
Umklapp process is relevant within the first or the third 
band, the aligned levels of the bands should favor the 
pairing, so that if the pairing correlation is suppressed, 
we can infer that the suppression is not an artifact. On 
the other hand, when the level of the second band is 
misaligned from Ep, one may naively think that the ef- 
fect of the Umklapp process within the band becomes 
obscured. However, the effect of the Umklapp processes 
should in fact be enhanced when the highest occupied 
level derivate from Ep, since the Umklapp processes be- 
come well-defined if the highest occupied level is doubly 
occupied. 



In the beginning we look at the half-filling (Fig. 17). In- 
deed, no enhancement of the pairing correlation is found 
and the over all decaying form is similar to the U = 
result as in the two-leg case at half-filling. 
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FIG. 17. Pairing correlation P{r) (□) against r for a 38 
rungs system for U = 2 with t± = 0.955 and 114 electrons 
(half- filling). The dashed line represents the noninteracting 
result. 



When n is decreased to make the second band half- 
filled, the enhanced pairing correlation is found (Fig. 18). 
Possibilities are either the Umklapp process is not rele- 
vant, or it is relevant but does not affect the other bands. 
In the latter density-wave correlation might be 

dominant due to a charge gap opening by the Umklapp 
scattering. However, at least in the sense of the weak- 
coupling theory, the charge gap in the second band only 
enhances the density-wave correlation at long distances 
from r~ 2 to r _1 (unity, the value of the exponent, comes 
from the gapless spin mode and it should be independent 
to U) in the weak-coupling limit and thus the pairing cor- 
relation may still remain dominant for small U. 
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FIG. 18. Pairing correlation P(r) (□) against r for a 38 
rungs system for U = 2 with t± — 0.87 and 110 electrons 
(the half- filled second band). The dashed line represents the 
noninteracting result. 

Although we cannot decide which of the above two 
possibilities applies, we do have a unique situation where 
the pairing correlation is enhanced despite the Umklapp 
processes being possible. This interesting situation does 
not appear in the two-leg ladder. 

Lastly, we study the case so that the first band is half- 
filled when n is further decreased and t ± is also decreased 
(Fig. 19). In this case, we again observe no enhancement 
in P(r), as expected from the weak-coupling theory and 
the above discussion. 
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FIG. 19. Pairing correlation P(r) (□) against r for a 38 
rungs system for U — 2 with t± — 0.725 and 74 electrons 
(the half- filled first band). The dashed line represents the 
noninteracting result. 

In this section, we have shown with the projector 
quantum Monte Carlo method that the enhancement of 



the pairing correlation expected from the results in the 
previous section is indeed found even for the interme- 
diate interaction strengths (U ~ 2t) and the interchain 
hoppings(ij^ ~ t). The features of the enhancement is 
similar to that in the two-leg case. 

We have also studied the cases where the Umklapp 
processes can be relevant. Especially, we found that the 
enhancement of the pairing correlation is not affected by 
the intraband Umklapp process within the second band 
which does not involve the pairing order parameter. 



V. PAIRING CORRELATION IN FINITE 
SYSTEMS — A CASE STUDY FOR THE ID 
ATTRACTIVE HUBBARD MODEL 

As we have seen in sections 2 (two-leg ladder) and 4 
(three- leg ladder), QMC result for the pairing correlation 
in Hubbard ladders are consistent with the weak-coupling 
prediction when the highest occupied level and the lowest 
unoccupied level (called the level offset Ae hereafter) in 
the relevant free-electron bands are made to lie close to 
each other. A similar situation is also found to occur in 
a QMC study of the single chain with distant transfers 
(t-t' -U model) H 

We believe that small level offsets should be required 
to mimic the thermodynamic limit (bulk systems) , where 
this quantity is infinitesimal. However, we have to quan- 
tify the criterion systematically. Otherwise, the results 
obtained in the previous section may be taken as an un- 
controlled finite size effect. In this section we actually 
give a circumstantial evidence for justifying that we can 
indeed quantify this, with which we can tell that the re- 
sults in the previous sections are meaningful. 

We start with asking ourselves the following question. 
Let us take a model (such as the attractive Hubbard 
model with an attractive interaction U < 0) that is 
known to superconduct. For a bulk system we know that 
the pairing correlation should have an slowly decaying 
asymptotic behavior for arbitrary [7^0, but in a nu- 
merical calculation such as QMC the tractable system 
size is limited, so that it is inconceivable that a slowly 
decaying asymptotic behavior is obtained even for, say, 
|J7| = O.Olt. So the question is what is the requirement 
for such an asymptotic behavior to be detectable in finite 
size systems. 

We shall conclude in this section that the level offset 
as compared with the spin gap is the key ingredient. The 
criterion will also shed light on numerical calculations for 
systems which may possibly have small spin gaps. 

As a case study, we take here the ID attractive Hub- 
bard model. This is bepause the model is exactly solvable 
with the Bethe Ansatzca, so that we know the exact form 
of the correlation function as well as the value of the spin 
gap for arbitrary values of parameters. We compare this 
with QMC results for the pairing correlation for system 
sizes exceeding one hundred sites. We deliberately choose 
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small attractions to look into the case where the spin gap 
is small. 

Let us first briefly recapitulate the exact result for 
the ID attractive Hubbard model, whose Hamiltonian 
is given by 

H = -t ^(c] CT Ci +lCT + h.c.) + u n n n n ( u < °) 

icr i 

(15) 

in standard notations. We consider the case where the 
numbers of up-spin and down-spin electrons coincide. In 
this case, the spin gap is present for all the values of 
U < and the band filling n. At half filling (n = 1), 
the pairing correlation decays as 1/r with the real space 
distance r, regardless of the value of U. Note that the 
power is reduced abruptly by unity from the power (=2) 
for U = as soon as an infinitesimal U is switched on. 
This is an effect of the spin gap, which opens as soon as 
an infinitesimal U is switched on, this is due to the fact 
that the spin phase is locked to give a long range order 
for its contribution to the pairing correlation function for 
U < 0. For a general value of n the pairing correlation 
decays as l/r 1 / Kp , where (1 <)K p (< 2), an exponent 
appearing in the Tomonaga-Luttinger theory, monotoni- 
cally increases with \U\, and decreases with n(< 1). 

We have adopted the projector Monte Carlo method to 
calculate the on-site singlet pairing correlation function 
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(16) 



in the ground state. The calculation is free from the neg- 
ative sign problem because we are considering the case 
of attractive Ul3 All the calculations are performed for 
a 114-site system. We consider the case where the closed 
shell condition (non-degenerate free-electron Fermi sea) 
is full-filled, so that there is a finite gap Ae between the 
highest occupied level and the lowest unoccupied level. 
To evaluate the decaying power of the pairing correla- 
tion, we least-square fit the numerical result to l/r a in 
the range r — 6 ~ 27 by taking logarithm of the data. 
In this range the periodic boundary condition is seen to 
have little effect. Although a cosine-like fluctuation is ob- 
viously present in the data, fitting all of the points in a 
finite range should average this effect out, namely, fitting 
cos(fcr)/r' 3 to a form a/r 1 would give a ~ 0. 

In Fig. 20, we first present results for half filling (n = 1), 
where the pairing correlation has an asymptotic behavior 
like ~ 1/r regardless ofU. Our QMC result for U = —t 
shows that the power is a ~ 1.4 in contradiction with the 
exact result. Only when U becomes as large as 2 do we 
recover a ~ 1. 
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FIG. 20. QMC result for the on-site pairing correlation 
function, P(r), plotted against the real space distance r in a 
114-site ID attractive Hubbard model at half filling. U — 
(dashed line), U = -t (O), U = -It (□). Straight dashed 
lines are least squares fit with l/r a for r — 6 ~ 27. 



Now, the spin gap Ag for this system size at U = 
—t(—2t) estimated from the Bethc Ansatz is 0.12t 
(0.23t). If we take the ratio of these values to the level 
offset Ae around the Fermi level, which is 0.1 It, they give 
values 1.1 and 2.1, respectively. This suggests that the 
pairing correlation, which is governed by the spin gap, 
cannot be clearly detected unless As/Ae is sufficiently 
large, in agreement with our speculation mentioned at 
the beginning of the present section. In fact Ag/Ae is 
infinite in a bulk system, and we also assume an infinite 
As/Ae in the Tomonaga-Luttinger theory as well. 

In Fig. 21, in order to confirm the point away from the 
half filling, we have also looked at n = 1/3, where As(= 
0.055£) is smaller than that for n = 1. The power for 
the pairing correlation is now a ~ 0.85 for U — —t form 
the Bethe Ansatz. The QMC result for U = —t gives 
a ~ 0.9, indicating that the numerical result is rather 
reliable already at U — — t in this case. If we look at the 
level offset, As/Ae ~ 1.6 is indeed greater than that for 
n = 1 and is closer to the case of U = —2t for n = 1. 
This result further confirms our view. 
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FIG. 21. QMC result for P{r) with U = -t. n = 1 (O), 
n = 1/3 (□). 



In conclusion, we have shown that the effect of the spin 
gap can be clearly seen in the correlation function of finite 
size systems only when the size of the gap is sufficiently 
large compared to the discreteness of the one-body en- 
ergy levels. Combining the results for the Hubbard lad- 
ders and for the t — t' — U model, we believe that this is 
the case in general: A criterion for reproducing the cor- 
rect form of the correlation functions is to make Ag/Ae 
large. This implies that to reproduce the correct asymp- 
totic behavior of the correlation functions in models hav- 
ing small spin gaps, very large systems are required. An 
alternative way is to tune, if possible, the values of the 
parameters so that Ae becomes small around the Fermi 
level, as we have done for ladders. 

One final point is, since Ae is defined as the level offset 
at U = 0, the quantity can become ill-defined for large 
interactions, where we have to consider that the quantity 
should be renormalized. 



VI. SUMMARY 

In the present paper, we have studied the pairing 
correlation in the Hubbard ladder with two(even) or 
three(odd) number of legs. This has been motivated from 
a conjecture due to Rice et al. that an even-numbered 
ladder should exhibit dominance of the interchain singlet 
pairing correlation as expected from the persistent spin 
gap away from half-filling. Naively, one can then expect 



that an odd-numbered ladder should not exhibit domi- 
nance of the pairing correlation reflecting the presence of 
gapless spin excitations. 

We have first considered the two-leg Hubbard ladder 
model. In the weak-coupling theory, in which the interac- 
tions are treated with the perturbative renormalization- 
group method, a d-wave like pairing correlation be- 
comes dominant reflecting a spin gap in the two-leg 
Hubbard ladder. The relevant scattering processes are 
the pair-tunneling process across the bonding and the 
anti-bonding bands and the backward-scattering process 
within each band. The pair-tunneling process is reminis- 
cent of the Suhl-Kondo mechanism for superconductiv- 
ity in the transition metals with a two-band structure. 
However, the weak-coupling theory is correct only for 
infinitesimally small interactions in principle. Thus the 
calculation for finite interaction U is needed, but existing 
numerical calculations for finite U have been controver- 
sial. 

In Sec. II, we have applied the projector Monte Carlo 
method to investigate the pairing correlation function in 
the ground state for finite U. We conclude that the dis- 
creteness of energy levels in finite systems affects the pair- 
ing correlation enormously, where the enhanced pairing 
correlation is indeed detected for intermediate interaction 
strengths if we tune the parameters so as to align the dis- 
crete energy levels of bonding and anti-bonding bands at 
the Fermi level in order to mimic the thermodynamic 
limit. The enhancement of the pairing correlation in the 
U = 2t case is smaller than that in the U = t case. This 
is consistent with the weak-coupling theory in which the 
pairing correlation decays as r~ x ^ 2K <^ [K p is the critical 
exponent of the gapless charge mode) at long distances 
and K p is a decreasing function of U . 

In the cases where interband or intraband Umklapp 
process is possible, the pairing correlation is not en- 
hanced. This result is also consistent with the weak- 
coupling theory. 

We then moved onto the correlation functions in the 
three-leg Hubbard ladder model in Sec. III. Whether the 
above 'even-odd' conjecture holds for the simplest-odd 
ladder (i.e. the three-leg ladder) with a plural number 
of charge and spin modes is an important problem. This 
has remained an open question, since there had been no 
results for the pairing correlation function in the three-leg 
t — J or Hubbard ladder models. 

A key is the coexistence of gapless and gapful spin exci- 
tations in the doped three-leg Hubbard ladder. This has 
been analytically shown from the correlation functions 
starting from the phase diagram obtained by Arrigonio, 
who enumerated the numbers of the gapless charge and 
spin modes with the perturbative renormalization-group 
technique in the weak-coupling limit. If we turn to the 
correlation functions, we have found that the coexisting 
gapful and gapless modes give rise to a peculiar situa- 
tion where a specific pairing across the central and edge 
chains (roughly a d-wave pairing) is dominant, while the 
2kp SDW on the edge chains simultaneously shows a sub- 
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dominant but still long-tailed (power-law) decay associ- 
ated with the gapless spin mode. The relevant scatter- 
ing processes are the pair-tunneling process between the 
top and the bottom bands and the backward-scattering 
process within the top and the bottom bands. The sit- 
uation is rather similar to the two leg case where the 
pair-tunneling processes play an impprtant role for the 
enhanced pairing correlation. Schulzc3 has independently 
obtained results for both the SDW and the pairing cor- 
relations which are similar to those given in Sec. III. 

However, as discussed in the two-leg case, it is not 
clear whether the weak-coupling results might be appli- 
cable only to infinitesimally small interaction strengths. 
In Sec. IV, we have thus checked the pairing correlation 
in the three-leg Hubbard ladder with the QMC method 
tuning the parameters so as to align the discrete energy 
levels of the first and third bands at the Fermi level as in 
the two-leg case. The enhanced pairing correlation is in- 
deed detected even for intermediate interaction strengths 
in the three- leg ladder. The enhancement of the pairing 
correlation in the U — 2t case is smaller than that in the 
U = t case as in the two-leg case in Sec. II. This result 
is consistent with the weak-coupling theory in Sec. Ill in 
a similar reason with that in the two-leg case. Namely, 
the exponent of the pairing correlation is a decreasing 
function of -Kp's which should decrease with U. 

Various effects of the Umklapp processes have also 
been discussed. Especially, it is found that the enhance- 
ment of the pairing correlation is not affected by the in- 
traband Umklapp process within the second band which 
does not involve the pairing order parameter. The key 
message obtained in the present paper, is that the dom- 
inance of the pairing correlation only requires the exis- 
tence of gap(s) in not all but some of the spin modes. 
This is independent to whether the number of legs is 
even(two) or odd(three). 

There are still important open questions. One is how 
the Hubbard model can possibly be related to real sys- 
tems such as the cuprates. Specifically, One should study 
the two-leg and/or three-leg Hubbard ladder for larger 
U, where the dominance of the pairing correlation might 
be lost. Furthermore, the three-leg t — J model should 
be studied to be compared with the Hubbard ladder. It 
would also be interesting to further investigate ladder sys- 
tems with larger numbers of legs. The Hubbard ladders 
with more large number of legs is also of interest, since it 
may have a key point to understand-,the two-dimensional 
systems. More recently, Lin et alr3 have examined the 
weak-coupling phase diagram for the four-leg Hubbard 
ladder where the pairing correlation is the most domi- 
nant in wide parameter region. 

Finally, in Sec.V, a QMC study for the pairing corre- 
lation in the ID attractive Hubbard model is given. The 
model can be exactly solved by Bethe Ansatz, so that 
we know the exact form of the correlation functions at 
long distances and the values of the spin gaps. A QMC 
study has shown that the effect of the spin gap can be 
clearly seen in the pairing correlation function of finite 



size systems and the behavior at long distances is consis- 
tent with the exact result only when the size of the gap 
is sufficiently large compared to the discreteness of the 
one-body energy levels. The result is a circumstantial 
evidence to justify the enhancement of the pairing cor- 
relation in Sec. II and Sec. IV and combining this result, 
we believe that this is the case in general for finite size 
systems. 

We are grateful to JProfeasor H.J. Schulz for directing 
our attention to refs.cl anoo. We wish to thank E. Ar- 
rigoni for sending us his work prior to publication. We 
also acknowledge helpful discussions with Professors M. 
Ogata, Y. Takada, H. Fukuyama, M. Imada, A. Fujimori, 
N. Mori, M. Azuma, Z. Hiroi, and J. Akimitsu. Nu- 
merical calculations were performed on a FACOM VPP 
500/40 at the Supercomputer Center, Institute for Solid 
State Physics, the University of Tokyo and on a HITAC 
S3800/280 at the Computer Center of the University of 
Tokyo. For the latter facility we thank Professor Y. 
Kanada for support in a Project for Vectorized Super 
Computing. This work was also supported in part by a 
Grant-in- Aid for Scientific Research from the Ministry of 
Education, Science, Sports, and Culture of Japan. One 
of the authors (T.K.) acknowledges the Japan Society for 
the Promotion of Science for a fellowship. 



APPENDIX A: CALCULATION METHOD IN 
SEC.III 

1. Derivation of eq. (3.4) 

Here we derive eq.(3.4) in the standard bosoniza- 
tion method. The intra- and the inter-band forward- 
scattering terms, which can be diagonally treated in the 
phase Hamiltonian as we will see in the following, are 
produced from the intrachain forward-scattering terms 
using eq.(3.3) as follows: 

2Trv F g x - x - t 

£ 2_^i 2—1 C fJ., + ,k 1 ,a C ^,-,k 2: cr' C tJ.,-,k2+cl,a' C ^, + M-q,<J 

ki ,k2,q /^,(T,er' 

= Hf + pair-tunneling terms, 
_ kv f g \ - v~^ro f t 

- 2 ^' ja i,+,fcl,ff a i,-,fc2,Q-' ai ,-,fc2+g,g ,0 l,+,fcl-g,<T 

ki ,/c2 : q cr,{T' 

+ 4a 2, + ,fc 1 ,o- a 2 -,fc 2 ,o-' a 2 ,-,fc 2 +9,(T' a 2, + ,fci- g ,cr 

I Q t t 

+ da 3,+,fc 1 , CT a 3,-,fe 2 ,er' a 3,-,fe2+9,o-' a 3,+,fei-g,ff 
+ 2 ( a l + ,fe 1 ,o' a 2 ,-,fc 2 ,<r' a 2 -M+q,a' a ^, + M-q,<? 
+ a 3, + ,fei,o- a 2 ,-,fc 2 ,cr' a 2 ,-,fe 2 +9,T' a 3, + ,fci-q,<T 
+ a L + ,k 1 ,<r a l-,k 2 ^' ai -M+q,<?' a 2, + M-q,<? 
+ a L + ,k 1 ,'r a l ,-,fc 2 ,cr' a 3 -,k 2 +q,a'a2,+,kx-<l,°) 
+ 3 ( a l,+,fc 1>CT °3 -,fc 2 , CT ' a 3,-,fc 2 +g )0 -'ai,+,fci-<Z,<7 
+ a 3, + ,fei,<r a l -,fc 2 ,cr' a l -M+q,"' a ^ + M-q,c)\ 
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+pair-tunneling terms. 



(Al) 



Here H{ consists of the forward-scattering processes in 
the band description and the electron operator with in- 
dex + (— ) belongs to the right- (left-) going branch. Wc 
prepare the following bosonic operators as in the usual 
single-chain case. 



/ \l/2 

(m ) a i.+, P -k,v a *,+,P,° for k > °> 

( « V'V ^ 

^ \k\L J l^i,p,cr <*»,- 



OCi,k 



,p+\k\,a 



IP,* 7 



for k < 0, 
(A2) 



1/2. 



— y 



(pfrj Ei, P , CT ^<-, p +| fe | ;(7 ai,-,P,a forfc<0. 

(A3) 



,t 



ai,+, P ,a for fe > 0, 



o-i,k(Pi,k) corresponds to the charge (spin)-dcnsity excita- 
tion in band i. Note that, a,^ and /^fc obey the boson 
commutation relation: 

[a,,fc,a:i',fc'] = [/3»,fe, /3»',fc'] = [a*,*;, A', fc'] = [a,,fc,/4,fc'] = 



a l',k'] = \Pi,k,0l t k>] = S i,i' 5 k,k>- (A4) 
iff is expressed in terms of only a, ; fe as 

Ht = ^rJ2 k [ 3 ( a \,k a i,-k + 4,fe4,-fe 

fc>0 

+ ai.fcai.-fc + a 3 ,fea3,-fc) + 4 04,k a 2,-fc + a 2,fcC*2 _ fe ) 

+ 2 ( a l,k a 2,-k + a i,fc a 2,-fc + a 3,fc a 2,-fc + a 3 : fc a 2,-fc 

+ a l,-fe a 2,fe + a l,-fc Q; 2,fe + a 3-k a 2,k + a 3,-A;a2,fc) 

+ 3 ( a l,k a 3,-fc + a l,k a 3,-fc + a3,fc a l,-fc + a 3,feQ!l,-fe)]- 

(A5) 

Furthermore we can rewrite the non-interacting part of 
the Hamiltonian Ho as 

Ho = v F J2 \p\ a l P a i,P + v fJ2 W\, P k,v ( A6 ) 

i,p i,p 

Now we introduce the phase variables as in the single 
chain case by the following equations: 



k>0 



-Afe/2 



ikx ( „ t 



x[e- lfe3: (aT fe ± aij _ fc ) - e lfa (a lifc ± a\_ k )\ (A7) 



-Afe/2 



fc>0 



x [ e -^-(/?l :fc ± ft,_k) - e^(A, fc ± ^_ fe )]. (A8) 



Here the phase variable 4>i+(9i+) can be regarded as the 
phase of the spin(charge)-density wave, while 0j_(#j_) is 
the field dual to 0j+(#j+). 

From (B.7) and (B.8), we can write the diagonal part 
of the Hamiltonian, Hd — Ho + Hi, in terms of the phase 
variables as 



& — Hspin ~t~ -^charges 

H spin = Y,^J <*s[(0*&+) 2 + {dx<k-f\, (A9) 

i 

^charge 

= £/ ^[(i+^)(^i+) 2 +(i-^)(^i-) 2 " 

+ S / da; [( 1 + 5)( 9 ^2 + ) 2 + (l-.g)(9^ 2 _) 2 " 

+£ / ^[( 1+ ^) ( ^ 3+)2 

(A10) 



+ 



^F5 
47T 



8tt 



+ (i-r ff ) ( ^3_)2 



dx[(a^ 1+ )(a^ 2+ ) - (0 x 0i_)(0 x 2 -) 

+ (5^ 3 +)(a^ 2+ ) - (a^ 3 -)(a^2-) 

da;[(a !B fl 1+ )(a x fl3+) - (OA-X^s-) 



Thus iJd is separated to both the spin-part H sv - ln and the 
charge-part I/charge • -^charge is also diagonalized by using 
eq.(3.5), while -ff sp in is already diagonalized. As a result, 
eq.(3.4) is easily obtained. 



2. Calculation of Correlation Functions 

Here we explain the method to derive the correlation 
functions. As examples, we here calculate the correla- 
tion functions of the intrachain singlet pairing in the edge 
chains and of the singlet pairing across the central and 
the edge chains. As stated in the text, the relevant scat- 
tering processes, are the pair-tunneling process between 
the first and the third bands and the backward-scattering 
process within the first or the third band. 

The pair-tunneling process is expressed in terms of the 
phase variables as follows: 

^3+1^3-1^1-4^1 + 1 + ^1+1^1-^1-1^1+1 

+ V4+^3-T^i+T^i--l + ^3+T^3-i^i+l^i-T + h - c - 
oc r]3+TV3-iVi-im+i;exp[i(8i- - #3- + <t>i+ - 03+)] 
+ %+i??3-T ? 7i-T ? 7i+iexp[i(6li_ - 6> 3 _ - </> 1+ + cf> 3+ )] (All) 

+ V3+lV3-lVi+lVi-ie x P[i(®i- - 03- + 0i+ + 03+)] 

+ Vd+^ilz-iili+iili-^vWi- - &3- - 0i+ - 03+)] + h.c, 

= 2cxp[z(6>i_ - 6> 3 _)][cos(0i+ - 3+ ) - cos(0i + + 3 +)] + h.c, 

= 8cos(V2xi-)sin(0i+)sin(0 3+ ). 
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Here we have defined the product of the U operators^ as 
Vi+^Vi-l — Vi+l 7 li-h but this convention does |»ot,affect 
the correlation functions as in the two-leg caseEaO. 

The backward-scattering process within band i is also 
expressed in terms of the phase variables through eq.(3.6) 
as 



Vi+lVi-lVi+lVi-iex-P 
-2cos(20;+). 



[-2i<j> i+ ] + h.c. 



(A12) 



This and eq.(All) give the eq.(3.7) in the text. In the 
beginning we calculate the correlation function of the in- 
trachain pairing in an edge chain (a chain). The order 
parameter is expressed in the band description as 

OintraSS = Ipa + ^a-l 

~ ^bPi+^i-l + ^3+1^3-1 + 2^2+T^2-l], (A13) 

where we have picked up only the two-particle opera- 
tors whose correlations show power-law decay at long 
distances. The correlations of the other two-particle op- 
erators correlation decay exponentially due to the gapful 
field(s). We can rewrite the above equation in the phase 
variables as 

^l+T^i-l + 2^2+T^2-j, + ^3+1^3-4, 

oc 77i +T 77i_;exp[j(6li_ + c/> 1+ )] + J7 3+T 773_j.exp[i(03- + 3+ )] 

+ 277 2+T 772-|exp[i(6l2- + (f>2+)], 

f... 1 1 1 , , o 

= icxp^iK— xi- + — X2- + —X3- 



V3 / 



V5' 



l l 



1 



X3- 



+-x P {i[(- ^ , v _ / 
+2iexp[i(0 2 - +0 2 +)], 

r. r/ 7T 1 1 , 7T-1 

f. r/ 7T 1 1 . 

+^H-2 + ^ X2 - + 7! X3 - )+ 2 1 l 

+2iexp[i(9 2 - - 02+)], 
= 2iexp[z(6» 2 _ - 02+)]. 



Now we calculate the correlation function of the singlet 
pairing across the central and the edge chains. The order 
parameter is expressed as 

OcESS = Oa+T + Ipl+lWp-l - Oa+1 + V+-+1 ) V>/3-T , 
~ Vl + T^l-l - -03+T"03-l 

-(^l+lV'l-T -^3+1^3-1)- ( Al6 ) 

Here again we pick up only the two-particle operators 
whose correlations show power-law decay. In terms of the 
phase variables, the order parameter can be rewritten as 

^l+T^i-J, - V's+T^-i - V'l+iV'i-T + V'3+i.V'3-T 
oc T}i + fqi-iexp[i(&i- + 0i+)] - 7? 3+T ?7 3 _jexp[i(03_ + 3 +)] 

- Vi+lVi-tex.p[i(^i- - 0i+)] + V3+lV3--[exp[i(6 3 _ - 3+ )], 
-ex P[i{ (-i=^ + i=X2- + ^X3-) + f}] (A17) 

r . r/ 1 7T 1 1 x TT-n 

- eX#{( Vl^ + ^ X2 -^ X3 - ) "2>] 

+exp[i{( -^^ + 7! X2 - + ^ X3 - ) -? }] ' 

r / 1 1 

= -4icxp[«(-^=X2- + -^X3-)J- 

Calculation of the interchain pairing correlation function 
is quite similar to that of the intrachain pairing correla- 
tion. 

(O C ess(z)Ocess(0)) 

( ex Pf*(^X2-(aO + -7=X3-( X ))] 



(A14) 



Here we have fixed 0i+ 



U3+ 



7r/2 and Xi- = ^/V^ 



as discussed in the text and the terms containing t he 
gapful fields are canceled out. Now we can calculate the 
correlation function: 



(O in traSs(a;)OintraSs(0)) 

oc (exp[i(6> 2 -(x) - 2 +(a:))]exp[i(02-(O) 
r 1 
2 

2tt 



exp[ 



2+(0))]), 

{<(&-(*) - £ 2 -(0)) 2 ) + ((02+(x) - 02+(O)) 2 }}], 



exp[ 



1 



3L^K* p2 



-(■ 



A' 



p3 



3)E : 



-Afc 



eXPh 6 ( A^ + A^ + 3)1 ° g(1+ A^ )] ' 

■T P 2 P 3 



vr 



>/6' 



xexp[i(i=X2-(0) + 4=X3-(0))]), 



1,1 



= exp[-^{^<(x2-(x)~X2-(0)) 2 ) 
+ ^((X3^W-X340)) 2 )}], 
X P 2 P 3 . 



(A18) 



From above calculations, we can see that the interchain 
pairing exploits the charge gap and the spin gaps to re- 
duce the exponent of the correlation function, in contrast 
to the intrachain pairing. 



(1 - cosfcc)], (A15) 
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